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Using a kinetic simulation of magnetic reconnection it was recently shown that magnetic-field- 
aligned electric fields (Ay) can be present over large spatial scales in reconnection exhausts. The 
largest values of Ay are observed within double layers. The existence of double layers in the Earth’s 
magnetoshere is well documented. In our simulation their formation is triggered by large parallel 
streaming of electrons into the reconnection region. These parallel electron fluxes are required 
for maintaining quasi-neutrality of the reconnection region and increase with decreasing values of 
the normalized electron pressure upstream of the reconnection region, /3eoo = S/roUeooTeoo/B^. A 
threshold (/3eoo < 0.02) is derived for strong double layers to develop. We also document how 
the electron confinement, provided in part by the structure in Ey, allows sustained energization 
by perpendicular electric fields {E±)- The energization is a consequence of the confined electrons’ 
chaotic orbital motion that includes drifts aligned with the reconnection electric field. The level of 
energization is proportional to the initial particle energy and therefore is enhanced by the initial 
energy boost of the acceleration potential, e$y = e Ayd/, acquired by electrons entering the 
region. The mechanism is effective in an extended region of the reconnection exhaust allowing 
for the generation of superthermal electrons in reconnection scenarios, including those with only 
a single x-line. An expression for the phase-space distribution of the superthermal electrons is 
derived, providing an accurate match to the kinetic simulation results. The numerical and analytical 
results agree with detailed spacecraft observations recorded during reconnection events in the Earth’s 
magnetotail. 


I. INTRODUCTION 


It is widely believed that magnetic reconnection ener¬ 
gizes electrons and is the source of superthermal electrons 
observed during solar flare events [I[ and through in situ 
measurements in the Earth’s magnetosphere [^. Several 
previous studies of reconnection conclude that the dissi¬ 
pation sites (where electrons are energized) are limited 
to electron scale regions and are insufficient in size to 
account for the large scale electron energization. To cir¬ 
cumvent this conundrum recent models have invoked the 
possibility of multiple reconnection locations and flux- 
rope formation, which can facilitate energy release in a 
larger spatial volume and thereby help energize the large 
number of electrons required by observations [^-Q. 

While there exists clear evidence of magnetic flux-ropes 
forming in the Earth’s magnetosphere, these flux-ropes 
do not appear to be volume filling such that energization 
at multiple x-lines becomes effective 0,0 . This suggests 
that energization within the exhaust from a single x-line 
needs to be efficient to account for the extensive energiza¬ 
tion levels recorded experimentally. We therefore revisit 
electron energization in single x-line reconnection and 
find that superthermal electrons are produced over large 
spatial scales in the single x-line reconnection exhaust 
when the total pressure of the plasma upstream of the 


reconnection region is dominated by the magnetic field 
iPeoo = “^UoPeoo/ ^1). In oui' kinetic simulation, a 
key to the strong electron energization is the break-down 
of adiabatic electron dynamics along the magnetic field 
lines, which, we find, is likely to occur when /3eoo < 0.02. 


The documented heating mechanism and non- 
adiabatic parallel electron dynamics lead to the forma¬ 
tion of strong magnetic field, aligned (parallel) electric 
fields, Ay. These structures are consistent with double 
layers, where local charge separation generates strong 
parallel electric fields Um- The double layers provide 
an initial energization of electrons streaming into the re¬ 
connection region (and exhaust) along field lines. In ad¬ 
dition, the double layer Ey helps confine the electrons 
within the reconnection region, allowing for further en¬ 
ergization by the reconnection electric field over the du¬ 
ration of multiple electron bounce orbits. The resulting 
electron distribution functions have characteristic signa¬ 
tures shaped by the interplay between the parallel and 
perpendicular electric fields during the energization pro¬ 
cess, which compare favorably to electron distributions 
observed in situ by spacecraft in the Earth’s magnetotail 
[3, [T3 - [l^ . As is the case for the magnetotail, /3eoo is also 
believed to be small upstream of reconnection sites asso¬ 
ciated with solar flares [l^, where the electron confine¬ 
ment provided by Eu may also be important to explain 
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the long lifetime of hard X-ray sources at flare loop-tops 

m r-i 

In previous work [IJ, [ij, have mostly an¬ 

alyzed scenarios including a guide magnetic field suf¬ 
ficient in strength that the magnetic moments, = 
mev']_l{2B), of the electrons are conserved, yielding adi¬ 
abatic (time reversible) electron dynamics for both the 
parallel and perpendicular electron dynamics. Here we 
consider anti-parallel reconnection, which includes un¬ 
magnetized chaotic electron behavior in the reconnection 
exhaust leading to irreversible electron dynamics. More 
important for the energization process, however, is the 
aforementioned break-down of the adiabatic electron dy¬ 
namics parallel to the magnetic field, causing the for¬ 
mation of large scale regions. The spatial extent of 
these regions is not tied to kinetic scales and expands 
with the reconnection exhaust. This suggests that the 
presented heating mechanism may be applicable to sys¬ 
tems like solar flares, which are much larger than the 
simulation domain [14| . 

The paper is organized as follows: In Sec. II we dis¬ 
cuss the formation of strong double layers, which form 
along the reconnection separators. The double layers are 
responsible in part for the confinement of electrons in 
the exhaust and we derive a threshold for their develop¬ 
ment. In Sec III we consider properties of the exhaust 
where pitch angle mixing is effective, and derive the rate 
for the energization of electrons trapped in the exhaust. 
In Sec. IV we investigate the electron distributions that 
form within the exhaust, providing an explanation for 
the observed electron flat-top distributions. In addition, 
a model is derived for the superthermal electrons which 
is in good agreement with the kinetic simulations results. 
In Sec. V the numerical and theoretical results are com¬ 
pared to spacecraft observations and the paper is con¬ 
cluded in Sec. VI. 


II. ELECTRON HOLES, AND DOUBLE LAYERS 

A. Electron pressure anisotropy, the acceleration 
potential, <1>||, and a condition for adiabatic electrons 

In this paper we document and analyze numerical 
simulation results obtained with the VPIC fully kinetic 
plasma simulation code [l^, implementing a configura¬ 
tion with /3eoo = 0.003. The two-dimensional simula¬ 
tion is initialized with a Harris neutral sheet, and open 
boundary conditions are employed in the x and z direc¬ 
tions. The initial particle distributions include counter¬ 
drifting Maxwellian ion and electron populations local¬ 
ized to support the Harris current layer. In addition, a 
separate background density is included, resulting in a 
total density profile n(z) = nosech!^(z/di) + Uoa, where 
no is the central Harris density, rioo is a uniform back¬ 
ground density and noo/nQ = 0.05 for this simulation. 
The ion-to-electron temperature ratio of the Harris pop¬ 
ulation is Tio/Teo = 5, whereas the temperatures, Tiao 


and Teoo, for the uniform background populations have 
the same ratio, but are a factor of three colder for both 
species. Lengths are normalized by the ion di = cjupi 
based on the Harris sheet density no, in a domain size 
of Lx X Lz = 320di X SOd^. Other parameters are 
milrfie = 400, T^oo = TOeC^/29 and Teoo = meC^/144, 
and ujpelLlce = 2, where Hce = ei?oo/(niec). The simula¬ 
tion data presented in Ref. were also from this run, 
but here we provide a more in depth analysis and refined 
understanding of the numerical results. 

Fig. [T] shows the profiles of the electron pressure 
anisotropy p||/pi and the acceleration potential for a rel¬ 
atively early stage of the reconnection process, t Ltd = 31. 
Here 4)|| is defined as 


$11 = 


E\\dl 


( 1 ) 


and is a pseudo potential which measures the integrated 
parallel electric field i?|| along the magnetic field lines. 
Because of the large electron thermal speed the magnetic 
field lines can often be considered stationary during a 
single electron transit [l^. Thus, e$|| characterizes the 
energy that electrons acquire in their free streaming along 
magnetic field lines into the reconnection region [20|. 

In Fig. [U the relatively large pressure anisotropy and 
large values of e^n/Teoo in the inflow are accurately de¬ 
scribed by the self-consistent adiabatic electron model 
derived in Refs. [ilHlil- The model provides equa¬ 
tions of state (EOS) for the electron pressure components 
P\\ = 4 * 11 ( 71 , 5 ) and pj_ = p±(ri,B). For anti-parallel re¬ 
connection the EOS are applicable to the inflow region 
where the electron magnetic moments, p = mev\/{2B), 
are conserved. The model includes the non-linear effects 
of electron trapping by the magnetic mirror force and by 
$ 11 . For the present geometry with e$||/reoo — 18, the 
trapped electrons dominate the properties of the electron 
fluid, and in this limit the electron pressure components 
P|| PC V? jB"^ and pj_ oc nB resemble the CGL scalings 
[2l|. In addition to the EOS, the model also yields the 
self-consistent values of $|| = $||(n, 5). 

Although the EOS do not apply to the electron dif¬ 
fusion region where the electrons become unmagnetized, 
through a momentum balance analysis of the electron 
layer it is shown that the EOS regulate the integrated 
current across the layer . Another result of the 

momentum balance analysis in Ref. ( 2 ^ is a scaling law 
for the magnitude of $|| expected just upstream of the 
electron diffusion region: 
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where the predicted values of e$||/Teoo become large 
at small values of /3eoo- Here h = n/n^ and jd^oo = 
‘2‘PonooTeoo /is the normalized electron pressure up¬ 
stream of the reconnecting current sheet. 

Beside the break-down of the magnetic moment as an 
invariant, a second mechanism also causes the adiabatic 
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FIG. 1: (a-c): Color contours of constant pressure anisotropy 
logic(Pll/Px)) s-od acceleration potential et&n/Teoo. The mag¬ 
netic field lines inside the island do not reach the ambient 
plasma and the acceleration potential, <E>||, is thus undefined 
within the areas left “white” in (b). 


EOS to become invalid. This second mechanism is re¬ 
lated to non-adiabatic effects in the parallel particle mo¬ 
tion, and occurs at small values of /3eoo ■ To estimate the 
critical value of /3eoo marking the transition to the non- 
adiabatic regime, we note that one important require¬ 
ment for adiabatic parallel behavior is that changes in 
$|j for a flux-tube moving into the reconnection region 
must be small during an electron transit through the re¬ 
gion. To quantify this condition we revisit the derivation 
of the adiabatic model in Ref. [l^, where a key element 
for solving the drift kinetic equation is an assumed or¬ 
dering 
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where we have used va — vjjL/D. For the low val¬ 
ues of j5eoo and assuming h = 1, Eq. m is approxi¬ 
mately e$||/Teoo = l/V/^eoo !)■ Furthermore, be¬ 
cause = rrie/(rriiPecx^), Eq. ([3]) may be written 

as /3eoo > \/TOe/wi, as required to ensure parallel adia¬ 
batic behavior. Thus, we expect non-adiabatic parallel 
behavior in the inflow (and along the separators) for 


(4) 

At the full mass ratio mi/me = 1836 this condition 
is j3eoo < 0.02, whereas for milm^ = 400 (applied in 
our simulations) the derived threshold is /3eoo < 0.05. 
Again, the simulation studied here is for a value /3eoo = 
0.003, and indeed, at later times in the run strong non- 
adiabatic behavior is observed, resulting in values of py 
and e<i)||/Teoo much larger than those predicted by the 
adiabatic theory. 



B. Electron holes, double layers and the formation 
of a large amplitude acceleration potential 
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Here D and L are the typical length scales across and 
along the reconnection region, respectively, with D/L ~ 
1/10. The ratio between the electron drift speed vp ~ 
O.luA and the thermal speed is small vp/vteoo — 1/400 
for P eoo — 1; but iucrcuscs with- dscrcusing vulucs of j^eoo- 

The adiabatic solution in Ref. [l^ corresponds to the 
limit where only electrons with small parallel energy may 
become trapped. This requires that during an electron 
transit time r ~ L/vtecx, the changes in the magnetic 
and electric well (trapping electrons) must be small com¬ 
pared to Teoo • If this condition is not satisfied the parallel 
electron behavior may become non-adiabatic. Following 
Eq. 13 of Ref. [l^, the condition for adiabatic dynamics 


As mentioned above, the acceleration potential with 
e4>||/reoo — 18 in Fig. [T](b) recorded for = 31 is 
consistent with the scaling law for <i>|| in Eq. ([2]). In con¬ 
trast, the acceleration potential for t^ci = 56 displayed in 
Fig. [3c) has a much larger magnitude, e^\\/Teos — 100, 
and is caused by the formation of structures in -Ey for 
which the largest amplitudes are observed within density 
cavities. The associated jumps in <1>|| are consistent with 
double layers iEl. The existence of such electron holes 
and double layers is well documented by spacecraft ob¬ 
servations within the magnetosphere. In the early obser¬ 
vations they were termed broadband electrostatic noise 
but Geotail observations showed they are solitons 
and not monochromatic mixture or coherent broadband 
tones [ 2 ^. More modern and detailed observations in¬ 
clude events recorded by the THEMIS mission [ 2 ^ 
and the Van Allen probes [ 2 ^. 

Two density cavities (with strong double layers) are 
clearly seen near the separators in Fig. H a) for x/di ~ 
200. In Fig. He) four points, d - g, are selected on a 
field line sampling the center of a large amplitude struc¬ 
ture in $ 11 . The corresponding distributions are shown 
in Figs. I3d“g)- For point d the value of e<i)y/Teoo — 5 
is relatively small; the distribution mainly consists of a 
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FIG. 2: (a-c): Color con¬ 
tours of constant den¬ 
sity logjQ(ne), pressure 
anisotropy logio(p||/pi), 
and acceleration potential 
e^ii/Te. Spatial loca¬ 
tions are marked in (c) for 
which /(u||,ux) is shown 
in panels (d) to (k). The 
magenta lines in (d) to (k) 
mark the trapped-passing 
boundaries as defined by 
Eq. (fT51l . 


beam of incoming electrons. Likewise for points e and 
/, the distributions are almost purely composed of the 
incoming beams energized by $|| with velocities reaching 
7 ri||/c ~ —1.2. For point g, in addition to the incoming 
beam, a population of electrons are observed mostly with 
velocities yu/c < 0.7. These relatively hot electrons orig¬ 
inate from the reconnection region, but given the large 
amplitude of $||, they are deeply trapped and cannot es¬ 
cape the reconnection region by parallel streaming along 
field lines. 

To provide more details on these large scale structures 
yielding direct acceleration of the incoming electrons, in 
Fig. [3l^a) a zoom-in-view is given for the region outlined 
in Fig. HJa). Within this region seven field lines are se¬ 
lected for which /||(x, U||) = 27 r f v±fdv± are computed. 
The resulting /|j(x, U||) are displayed in Figs. ISJb-h) as 
functions of x along the selected field lines. The upper 
most field line (b) in Fig.[31[a) is just upstream of the den¬ 
sity cavities. The corresponding /|| along this field line is 
shown in Fig.[31Jb). Signatures of developing instabilities 
are seen for x/di ~ 170. 

Strong instabilities are visible along the second field 
line (c) for which /|| displays evidence of both electron 
holes and double layers. The electrons with yun > 0 are 
streaming away from the reconnection region. The first 
double layer for x/di — 170 is seen to help confine elec¬ 


trons energized in the X-line region. However, the second 
double layer at x/di — 200 develops without the presence 
of energized electrons from the reconnection site. 

Figs. [31[d,e) illustrate how the double layers continue 
to develop. The vertical magenta lines marked d - g 
in Fig. [3je) correspond to the distributions in Figs. [UJd- 
g), respectively. The beams in Figs. [IJd-g) are readily 
identified as the structure (with < 0) in /|| containing 
incoming electrons continuously accelerated by $||. For 
these held lines at the center of the density cavity, e<i)|| is 
larger than the energies of electrons streaming away from 
the reconnection site. Thus, for x/di > 200 almost no 
escaping electrons (electrons with ynn > 0) are observed. 

The magnitude of the double layer in Fig. |3Kf) is still 
large. In addition to the incoming beam, the double 
layer is now being hlled with energetic electrons stream¬ 
ing away from the reconnection site and rehected back 
again by the double layer. Considering the evolution of 
/|| in Figs.[3Kb-f) it is clear that the acceleration by $|| of 
the incoming electron beam and the subsequent mixing 
in electron holes is responsible for the main heating of 
the electrons on these held lines intersecting the density 
cavities. 

Deeper into the exhaust the value of <i)|| is reduced. 
This is evident in Figs. [3Kg,h) where the incoming elec¬ 
tron beam now has velocities on the order of 7 n||/c ~ 
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FIG. 3: a) Color contours of constant 
density, log^Q(ne). Strong density cav¬ 
ities are observed where the density is 
reduce by a factor of 10 below the up¬ 
stream value. Seven field lines are se¬ 
lected and marked h to h. Along each 
field line the parallel distribution is com¬ 
puted, /||(u||) = 27r J fv±dv±, and dis¬ 
played in panels b) to h) as functions of 
7 U|| and x/di. The magenta lines in e) 
marked d-g indicate the positions corre¬ 
sponding to the distributions in Fig.[2jd- 
f), respectively. 


—0.8. The reduced values of $|| allow enetgized electrons 
to escape the reconnection region and finite values of f\\ 
are thus observed for x/di > 200 with 7 U|| > 0. These 
escaping electrons carry a significant heat flux away from 
the reconnection region. 


C. Formation of strong double layers 

Recently, Li et al. investigated the forma¬ 

tion of double layers at the interface between hot elec¬ 
trons energized at the reconnection site and cold elec¬ 
trons streaming into the reconnection region along field 
lines. Their numerical studies document electron holes 
and double layers with signatures similar to those ob¬ 
served in our simulations outside the density cavities 
(i.e. in Figs. [3ljg,h)). They modelled the hot electrons 
from the reconnection region as a Maxwellian with a 
temperature Th and observed double layer amplitudes, 
~ 0.73T;i, less than Th- This moderate amplitude 


allows for a significant fraction of the energetic electrons 
to escape 0 - 

Interestingly, the strong double layers that form within 
the density cavities in Figs. ISJc-f) appear to have some 
characteristics different from those studied by Li et al.. 
For example, we here observe an amplitude much larger 
than the energy of electrons energized at the reconnection 
site, allowing no electrons to escape by parallel stream¬ 
ing along field lines. Furthermore, the development of 
these large amplitude double layers is not driven by hot 
electrons streaming away from the reconnection region. 
This is clear because the strong double layer structure 
at x/di — 200 in Figs. ISjc-d)) develops before energetic 
electrons from the reconnection region reach this loca¬ 
tion (about 40 di downstream of the x-line). Below we 
argue that the strong double layers (distinct from the 
double layers outside the density cavities) form not just 
to reduce free streaming losses, but primarily to boost 
the electron density within the reconnection region. 

From the principle of quasi-neutrality, electric fields in 
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a plasma develop to maintain near equal densities of elec¬ 
trons and ions. To elucidate the mechanisms driving the 
strong double layers, in Fig. SJa) contours of the paral¬ 
lel electron flux are shown for the reconnection region, 
and large fluxes of electrons \neUf.\\\/nooVteoo — 4 flow to¬ 
ward the reconnection region. Thus, this flow is stronger 
than the incoming parallel flux of upstream electrons 
Tiioo = 27r v±dv±v\ifoo = nooVteoo/ 

To understand why the parallel flux of electrons is 
needed, we consider the blue area, Ai{= Ldl) on the 
schematic flux-tube in Fig. IHb). During the reconnec¬ 
tion process this area convects into the larger green area 
A 2 . From magnetic flux conservation within the flux- 
tube, the areas of the blue and green regions are related 
as 



where Bi and B 2 are representative magnetic field 
strengths within the two areas. As the in-plane mag¬ 
netic field vanishes at the X-line, as is ubiquitous in re¬ 
connection we have Bi > i? 2 . The areas considered are 
within the ion diffusion region where the ions are decou¬ 
pled from the magnetic field lines. The simulation shows 
that within this region the ion density, n, is nearly uni¬ 
form, such that area A 2 includes an increased number of 
ions compared to area Ai. This increase can be estimated 
as 

AN = n{A 2 — Ai) = nLdl ~ • (5) 

Meanwhile, the electrons are frozen in to the magnetic 
field in their perpendicular motion, so the matching in¬ 
crease AN of electrons (required for quasi neutrality) 
must be supplied by a parallel electron flux neUe\\. 



FIG. 4: a) Contours of parallel electron flows into the re¬ 
connection region, b) Schematic illustration of a flux tube 
section (blue) expanding as it enters the reconnection region 
(green). To maintain a uniform electron density (required for 
quasi neutrality) electrons flow along field lines into the green 
region driving the formation of the density cavity marked in 
red. 


To estimate Ue\\ we consider again the characteristic 
length scales D and L across and along the reconnection 
region. Previous studies (including Ref. [3l| ) have estab¬ 
lished that the inflow velocity upstream of the ion dif¬ 
fusion region is {D/L)va, with DjL ~ 0.1. Meanwhile, 
inside the ion diffusion region the velocity of the flux- 
tubes are enhanced by the factor Bi / B 2 yielding z;in,e = 
vaBiD/{LB 2 ). Thus, region Ai with length L convects 
into region A 2 during a time r ~ D/vin,e = LB 2 I(vaBi). 
The increase in the number of electrons within the con¬ 
sidered area can then be estimated as AN = nUe\\dlT. 
Equating this with Eq. the factors of n, L and dl 
cancel such that 

Bi f Bi \ 

Ue\\ = VA^ - Ij ~ dUA , 

where we have assumed a ratio — 2.5. This result 

may be rewritten on the form 



so for the present simulation (and in agreement with 
Fig.HKa)) we then obtain Uf,\\ — ‘ivteoo- 

In kinetic and also Hall MHD simulations of recon¬ 
nection at larger values of Peoo less pronounced density 
cavities are typical along the separators, and it has been 
argued that they develop to maintain pressure balance 
perpendicular to separators [ 3 ^ . We here emphasize the 
role of electrons streaming into the reconnection region 
as the more important cause for the near depletion of 
electrons from the density cavities. To the best of our 
knowledge, this additional mechanism for density deple¬ 
tion is described here for the first time. We also note 
that the resulting pattern of the parallel electron cur¬ 
rents along the separators are largely responsible for the 
characteristic Hall magnetic field of the reconnection re¬ 
gion [13 ■ 

The development of the observed strong double lay¬ 
ers is likely to be suppressed when nooUe\\ < Fnoo = 
«-ooi’teoo/(2i/7r)j allowing thermal streaming of electrons 
to supply the reconnection region with the electrons 
needed. With Eq. ®, a crude criteria for the devel¬ 
opment of the strong double layers is then 

/3eoo< 200^ . (7) 

rUi 

For mi /me = 400 this criteria suggests that the strong 
double layers may already form at /3eoo <0.5. However, 
the derivation of Eq. © does not include the effects of 
the adiabatic $|| given in Eq. o helping to boost the 
electron flow into the reconnection region. In contrast, 
the condition in Eq. ([4]) includes the flows driven by <l?|| 
and therefore provides a more accurate threshold for the 
transition to the non-adiabatic parallel dynamics. Nev¬ 
ertheless, we include the above derivation of Eq. © as it 
elucidates the mechanisms we believe are important for 
the formation of the density cavities and the associated 
strong double layer formation. 
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FIG. 5: a,b): Contours of con¬ 
stant density for tQ.ci = 56 and 
tQ.ci = 61, respectively. c) 

Typical electron test orbits tra¬ 
jectories. In the left side ex¬ 
haust the electron magnetic mo¬ 
ments are conserved, whereas 
the chaotic electron motion in 
the right side exhaust is non- 
adiabatic. d,e): Electron tem¬ 
perature components evaluated 
along the lines shown in b). 


III. PROPERTIES OF THE PITCH ANGLE 
MIXED EXHAUST 

A. Signatures of the pitch angle mixed exhaust 


The dynamics in the right side exhaust (for x/di > 
160) in Fig. [5] are different from the left side not only 
because of the asymmetry in $||. Another main differ¬ 
ence is that within the right side exhaust the magnetic 
moments of the electrons are not conserved when the 
electrons cross the mid-plane. The distributions for the 
points marked /i - fc in Fig. He) are shown in Figs. Hh- 
k). The points h and i are both at the mid-plane and 
we notice how their distributions are independent of the 
pitch angle 6 — Z(v, B). This is the signature of com¬ 
plete pitch angle mixing. 

The magnetic moment, /r = mv\l2, is only an adi¬ 
abatic invariant of the electron motion when the ra¬ 
dius of curvature of the magnetic field Rc is larger than 
the electron Larmor radius pi = mu_L/(ei3). This re¬ 
quirement for adiabatic motion may also be expressed 
as = Rc/Pi > 1 along the full electron trajectory 
34| . In the center of the exhaust the strongly bent field 
lines in combination with the low magnetic field strength, 
can lead to non-adiabatic electron motion with < 1. 
This causes electrons to pitch angle mix, washing out 


anisotropic structures in velocity space. 

To further explore the temporal evolution of the 
plasma, in Figs. Hs-jb) density profiles for tVLd = 56 and 
tU,ci = 61 are given, respectively. The main density cav¬ 
ities of tSlci = 56 are seen to move downstream with 
the exhaust while new cavities form closer to the x-line. 
During this downstream propagation the cavities remain 
characterized by large values of , such that the region 
of parallel energization is expanding in time. 

Examples of test particle orbits are given in Fig. [S](c). 
For the part of the blue trajectory which falls within the 
left side exhaust, the magnetic moment is an adiabatic 
invariant. Meanwhile, for the red trajectory typical of 
the right side exhaust, the magnetic moment is not con¬ 
served where the particles pass through the midplane of 
the domain. 

The described particle motion has direct implications 
on the pressure profiles. This is evident in Figs. [SKd,e) 
where the parallel and perpendicular pressure compo¬ 
nents are evaluated along the lines marked d and e in 
Fig. Hb). Consistent with the EOS by Le et al. for both 
inflow regions we observed Tc\\ » Te±. At the end of 
the electron diffusion regions pitch angle mixing leads to 
Te|| = Te±. Eor the left side, the magnetic islands lead to 
an increase of the field strength in the exhaust such that 
the magnetic moments of the electrons are conserved, and 








































































the betatron heating is effective. This leads to the tem¬ 
perature components with Te_L > Te||, and is in contrast 
to the right side exhaust where pitch angle mixing yields 

Te± — Tell. 


B. Simulation profiles of the pitch angled mixed 
exhaust 



2IK) I’’ x/d^ 240 


FIG. 6: Kinetic simulation results showing contours of a) 
e<l>||/Te, b) Te = rr(Pe)/3ne, c)14|| , d) logio(B), and e) 
Usuper/n representing the density of electrons with yv/c > 
2.5. f) parallel phase-space distribution /||(7r>||) as a function 
of X along the field-line highlighted in white in a)-e). The 
simulation results shown correspond to time tflci = 61. 

Key profiles for the electron dynamics at = 61 are 
given in Fig. [51 The area of a significant acceleration po¬ 
tential (e<l>||/Teoo > 30) in Fig. [SDa) is much larger than 
at previous time, tilci = 56, considered in Fig. HDc). In 
the right side exhaust, the boundaries for large values of 
e$||/T'eoo (> 30) coincide with the boundaries where the 
effective electron temperature, Te = (p|| -I- 2p±)j3ne, in 
Fig. HDb), displays an abrupt increase (by about a fac¬ 
tor of 20) from its value in the inflow. This is consistent 
with our conclusions above that the of the electron 
holes and double layers provides strong energization of 


the incoming electrons. This parallel heating is most sig¬ 
nificant within the density cavities, which, as displayed 
in Fig.jSDc), are characterized by the strongest streaming 
of electrons directed towards the reconnection region. 

As discussed above, due to the low values of B observed 
in Fig. inKd) along the exhaust mid-plane, the electron 
motion in the right side exhaust is non-adiabatic. Al¬ 
though the pitch angle mixing is limited to this region 
of low B, due to the non-localized motion of the trapped 
electrons, the pitch angle mixing impacts the electron dis¬ 
tributions in the full width of the exhaust. We also note, 
that in contrast to the bulk electron heating, the genera¬ 
tion of superthermal electrons occurs gradually in the ex¬ 
haust. This is evident in Fig. me) where density contours 
are given for superthermal electrons with yu/c > 2.5. 

The parallel electron dynamics along a typical exhaust 
field line is illustrated in Fig. mO- in Fig. |3J the con¬ 
tours represent constant values of /||(u||) = 27r f fv±dv±, 
here as a function of x along the field line highlighted 
in white in Figs. Ha-e). Again, large electron hole 
structures are observed, especially at the interface (for 
x/di ~ 240) between the cold incoming electrons and 
electrons already heated in the exhaust. 


C. Relative importance of Ay and E± for electron 
energization in the pitch angle mixed exhaust 

It is interesting to explore the relative importance of 
parallel and perpendicular electric fields for energizing 
the exhaust electrons. The local energy gain of the elec¬ 
tron fluid is quantified by E • Jg. As shown in Fig. 
we divide a part of the simulation domain into sectors 
of equal increments in the in-plane magnetic flux. Ay. 
The energy exchange terms Ty Je||, E • Jg, and E • are 
integrated over each sector and displayed in Fig. |7Kb)- 
Within the sector including the separator and the fol¬ 
lowing two downstream sectors, Ty Je|| amounts to about 
50% of the total E • Jg. As shown in Fig. me); these sec¬ 
tors include regions of strong electron flows Egy toward 
the X-line. Because <i>|| is generally increasing as the X- 
line is approached, incoming electrons are energized by 
E||. 

Further downstream of the X-line the electron flow is 
away from the X-line (see Fig. HKc))) such that the elec¬ 
trons endure a net loss of energy through work against 
E||. In fact, the average value of Ty Jgy in the exhaust is 
small, while the heating by Ej^ • Je_L accounts for nearly 
all the electron energization when integrated over the ex¬ 
haust. Also noteworthy, as can be seen in the traces in 
Fig. mb) the electrons acquire about 20% of the total 
dissipated magnetic energy, which is large compared to 
levels typically observed in simulations at higher /3eoo- 
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FIG. 7: a) The cross-section is sectored into intervals of 
Ay. For each sector i5|| Je||, E • Je, and E • Ji, is evaluated as 
displayed in b). 


D. Rate of electron energization by the 
reconnection electric field 


It is natural to suspect the reconnection electric field, 
iilrec, to be the cause of the documented energization 
by Ej^. The rate of this energization is given by P = 
ErecJey, where Jey is the electron current in the y- 
direction. The strongly curved field lines of the exhaust 
are favorable for energization related to the curvature 
drift [ 3 ^. This is shown directly in appendix A where 
we apply the guiding center model and obtain a Fermi- 
like energization rate proportional to the particle energy. 
However, it is not obvious that the guiding center model 
is applicable to the present exhaust where the magnetic 
moments of the electrons are not conserved. This moti¬ 
vates a more general derivation of the heating rate with¬ 
out the assumption of adiabatic electron motion. There¬ 
fore, we will estimate the heating rate using a more gen¬ 
eral formulation. This estimate will also include the pos¬ 
sible contributions from the magnetization currents (nor¬ 
mally neglected in calculations based on the guiding cen¬ 
ter model [ 3 ^). The end result is the confirmation of the 
heating rates obtained with the guiding center model. 

Our generally valid approach for estimating the heat¬ 
ing rate of electrons in the pitch angle mixed exhaust 
is first to generalize the MHD-force balance equation, 
X B 5 = Vp 5 , where the subscript <5 denotes a par¬ 
ticular class of electrons with ~ 'uq for t = to- We 
may thereby determine Jsy as a function of the electron 
energy, which in turn will provide us with the energiza¬ 
tion rate due the to the reconnection electric field as a 
function of energy. 

To limit the analysis to electrons with an initial veloc¬ 
ity around vg at time tg, we introduce /^(x, v,t) which 


at t = to is given by 

/(x,v,t) , -uo - < |v| < uo + 

0 , elsewhere 


/5(x,v,t) = 


( 8 ) 

Naturally, the time evolution fg is governed by the Vlasov 
equation 

^+v-V,/5-— (E + vxB)-V,/5 = 0 , (9) 

ot rrie 

and moments over fs are defined in the regular fashion 
as 


{^‘)s = J ^^fsd^v , 

yielding ng, and for / = 0,1, and 2, respectively. 

In analogy with the derivation of the standard electron 
momentum equation, by taking the I = 1-moment (or the 
v-moment) of Eq. (jS]) we obtain a momentum equation 
for the selected electrons 

meUg^^ = -eng{E -|- x B) - Vpg - V -ng . 

at 

Here we have split the pressure tensor, P^, into its scalar 
part pgl and the shear stress irg. For the pitch angle 
mixed exhaust, the latter is assumed to be negligible, 
TTi ~ 0. For t = tg, all electrons contributing to pg have 
velocities ~ uq and we thus find pg = (2me/S)ngVg = 
{4:/3)ng6g, with £g = 

Our aim is next to derive an expression for the average 
electron drift in excess of the E x H-drift. We therefore 
introduce — E x B/H^. Neglecting inertia, 

TO —>■ 0, we obtain the following momentum equation for 
electrons within the considered velocity interval: 

4^ 

0 =-en5 (E||-I-ud5 X B) - ^^”-5 ■ 

The current density carried by these electrons is = 
—enguug and the momentum balance perpendicular to 
the magnetic field is then governed by 

4:£n 

Ji X B = —V±ng . (10) 

An approximate heating rate is readily obtained from 
the above equation. The average z component of the left 
hand side may be estimated as 

{3g X B), ~ {Jgy) {B) , (11) 

where (...) denotes an average over the z-direction. Fur¬ 
thermore, to estimate the right hand side of Eq. m 
we write Vng ~ ngjD such that Eq. (HOI) approximately 
yields 

(Jjy) ^ 4 £g 

ng 3D (B) 
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Next, with Spec — VA {B) we obtain the following result 
for the heating rate 

dS {Jsy} -^rec ^ 
dt ns 3D 

A comparison with the simulation shows that this ex¬ 
pression overpredicts the heating rate by about a factor 
of two. However, using the numerical profiles to improve 
the various estimations above (such as ~ nsjD and 
Eq. (fTTll l we obtain our final expression for the heating 
rate consistent with the simulation 


d£ ^ ^ VA 
dt ~ °2D 


( 12 ) 


Accordingly, electrons trapped in the exhaust double 
their energy during an Alfvenic transit time across the 
full width of the exhaust, 2D. 

The energization rate in Eq. (fT^ is proportional to the 
initial energy of the electrons, and is similar to the rate 
for electrons in contracting magnetic islands Q. This 
magnetic island model was introduced to allow for en¬ 
ergization in extended reconnection exhausts filled with 
a bath of magnetic islands. Meanwhile, according to 
Eq. m electron energization in the reconnection ex¬ 
haust is effective even without magnetic islands and is 
not contingent on the development of pressure anisotropy 
(the latter has been found essential for energization by 
magnetic island @). In our analysis, the only require¬ 
ment for electron energization is that electrons are con¬ 
fined to the reconnection exhaust, which here is aided by 
the strong $|| associated with the double layer formation. 

We further note that for a magnetized exhaust, our 
starting point in Eq. m can be derived directly by in¬ 
tegrating the current carried by the guiding center drifts 
vjy and magnetization current 3e± = e f vof d^v -I- V x 
M, with the magnetization eM = bf jif d^v. Thus, the 
present analysis is therefore consistent with similar find¬ 
ings obtained using the guiding center model [s^. In fact, 
in Appendix A we show that Eq. (HU is readily derived 
using the guiding center model. 


IV. FLAT-TOP DISTRIBUTIONS AND 
SUPERTHERMAL ELECTRON ENERGIZATION 

In this Section we explore the heating of electrons in 
the exhaust. A statistical model is derived which charac¬ 
terizes the energy distribution of the superthermal elec¬ 
trons. While these electrons are mainly energized by 
perpendicular electric fields, we document how the ini¬ 
tial energization and confinement provided by <I>|| signifi¬ 
cantly enhances the effectiveness of the electron energiza¬ 
tion process. 


A. The trapped passing boundaries 

In Figs. [5Ka,b) we consider different classes of electron 
trajectories reaching the points highligted. The trajecto¬ 


ries passing through the reconnection region along field 
lines without any reflections are marked by A and D, and 
we denote these as passing. Meanwhile the trajectories 
marked B and C we denoted as trapped, with electrons 
bouncing back and forth along field lines, as the field 
lines convect across the reconnection region. The four 
classes of trajectories divide the (u||, u_L)-plane as shown 
in Fig. [HKc). The trajectories of velocity regions A and 
B reach the point considered with a negative value of U||, 
whereas uy > 0 for trajectories of regions C and D. The 
boundary between regions A and B and the boundary 
between regions C and D we denote the trapped/passing 
boundaries. In in Fig. |5Kc) these are shown by the ma¬ 
genta lines. 

When the magnetic moment is conserved the perpen¬ 
dicular energy is given by = fiB, implying that the 
parallel energy is £\\ = £ — fiBao- It then follows that 
the trapped/passing boundaries are found by solving the 
equation 

£\\oo = £ - e^W - fJ-B^ = 0 , (13) 

which expresses the physical condition that marginally 
trapped electrons will deplete their parallel energy 
(f||oo = 0) as they barely escape along the magnetic 
field lines away from the reconnection region. We note 
that for ux = 0 the trapped/passing boundaries start at 
kill = Vrj) = ^2e^\lme and at large v they asymptote 
to the line in Fig. Elc) characterized by the angle a, with 
cos(a) = B/Boo- 


B. Size of the loss-cone in velocity space 


For the analysis below it is convenient to introduce 
Rioss as a measure of the size of the electron loss-cones, 
where 47ri?ioss is the solid angle in velocity space of the 
loss-cones. Using the expression for the trapped/passing 
boundary in Eq. (IT^ it is readily shown that 


Ri 


OSS 


= 1 - 


1 - 


B 

Boo 



1/2 


Introducing the characteristic velocity of the acceleration 
potential = yt2eWi^Jrn^, an approximate form is ob¬ 
tained for the limit H <C Uq 


valid for |v| > v^. For |v| < the confinement is abso¬ 
lute with i?ioss = 0. 


C. Generation of flattop distributions in the pitch 
angle mixed exhaust 

The bulk exhaust distributions for the time slice of 
Fig. El are displayed in Figs. El First, consider the distri¬ 
bution in Fig. ETa) corresponding to the point marked 
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a 



Trapped - Passing 



FIG. 8: a) Classification of electron trajectories reaching a 
point in the inflow region. Trajectories A and D are passing 
while trajectories B and C are trapped, b) Similarly to the 
inflow region, the trajectories in the exhaust can be classified 
as passing and trapped. In the exhaust of antiparallel recon¬ 
nection, pitch angle mixing always occurs at the end of the 
electron diffusion region. In some cases, pitch angle mixing 
remains effective far downstream of the reconnection site, c), 
Regions in the (u||, u^)-plane identifying the velocity regions 
of the four classes of electron trajectories. The magenta lines 
represent the trapped/passing boundaries of Eq. (1 1311 . 


“3” in the panels of Fig. [51 As above, the velocity 
space (spanned by yrin and is divided into regions 
A, B, C, and D with distinct properties of the associ¬ 
ated electron trajectories classified in Fig. [5] Again, the 
trapped/passing boundaries between regions A/B, and 
regions C/D, are computed using Eq. (fT^ . Trajecto¬ 
ries of velocity region A are incoming electrons stream¬ 
ing along field-lines toward the reconnection layer. The 
electrons in velocity region C and D are all streaming 
away from the X-line, but the electrons in region C are 
trapped and will reflect into region B. Meanwhile, the 
electrons in region D have sufficient parallel energy to 
escape and will exit the simulation domain at the open 
boundaries. 

The changes in the Fig. [HKa) distribution are dramatic 
across the A/B boundary calculated using Eq. (fT^ . The 
region A electrons are well characterized as beams with 
a parallel energy e$||. A small spread is observed about 
the center of the beams corresponding to the low up- 



FIG. 9: Golor contours of the electron distribution function 
for the points marked P3, P4 and P5 in Fig. [6] For the distri¬ 
butions at P3 and P4 the velocity regions AtoD correspond 
to the distinct trajectories shown in Fig. [TKa) and described 
in the text. Because of the low values of the magnetic field in 
the center of the exhaust, pitch angle mixing is effective and 
the P5 distribution is fully isotropic. For energies £ < e4>|| 
the values of / are nearly constant, resembling flattop distri¬ 
bution observed in situ in the magnetotail. 


stream temperature of these incoming electrons. The dis¬ 
tribution in Fig. |9l(b) is for the point below the midplane 
marked “4” in the panels of Fig. [51 Given the chances 
in direction of B^, we note how these beams of incoming 
electrons have opposite signs of yny above and below the 
mid-plane. 

The region A beams terminate at the midplane of the 
simulation domain where, as shown in Fig.lHKd), the mag¬ 
netic field is weak, B/Bq < 0.1. This is also evident in 
the distribution in Fig. ETc) corresponding to the point 
marked “5” in Fig. [6] on the same field-line as points 3 
and 4. The weak magnetic fields causes a complete break¬ 
down of the magnetic moment as an adiabatic invariant 
and the resulting chaotic electron motion causes the dis¬ 
tributions within this layer to become fully isotropic in 
velocity space. 

As the isotropized electrons of velocity region A leave 
the midplane they populate velocity regions B, C and 
D throughout the exhaust. In principle, these elec¬ 
trons should only add to / at their injection energy 
£ = e$||. However, distributions in velocity space with 
negative slope {df/dv < 0) are unstable and instabili¬ 
ties (including electrons holes) develop in the simulation 
which are responsible for scattering of the electron en- 
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ergies. These processes and the energization mechanism 
described by Eq. © drive the distribution function to¬ 
wards df/dv ~ 0, such that for £ < e$[| flat-top distri¬ 
butions are approached where / ~ const [l3,113, ■ 


D. Statistical model for superthermal electrons 

We now seek to develop a statistical model for the 
electron energization, which includes the free streaming 
losses for electrons at energies above e$||. Given the pitch 
angle mixing in the exhaust, the phase space density is 
assumed to be fully isotropic /(v) = f{v). To obtain an 
evolution equation for f{v) we consider the phase-space 
volume bounded by vi and V 2 in Fig. 1101 This spheri¬ 
cal shell in velocity space has a width dv = V 2 — vi and 
a volume Airvldv. The number of electrons in this vol¬ 
ume changes because of three effects. 1) Electrons within 
the losscone of size AirvidvRioss are all lost within half a 
bounce time Tb/2. 2) Electrons with velocities just below 
V 2 will be lost by acceleration across the v = V 2 bound¬ 
ary. 3) Similarly, electrons with velocities just below vi 
will be accelerated into the volume. Considering a time 
interval dt, particle conservations may then be expressed 
as 


Airv'^dv -;^dt 
dt 


= — A-Kv^dv f Rio 


dt 


’ n /‘2 


- 4^Atf 
dt 


+ dTT V^^dtf 
dt 


It is clear that the last two terms in the equations above 
can be rewritten as dv times a differential. Dividing 
through by 47r dv dt we then obtain. 


V 


2 


dt 


^fR 

J -^\o 

n 


dv 



(15) 


To proceed requires estimations of Tf, and dvjdt as 
functions of v. From the trajectory on the right hand 
side of Fig. [SKc) we conclude that the characteristic orbit 
length is about 2D, where (as above) D is the half width 
of the reconnection exhaust. Then, the orbit bounce 
time is approximately ~ 2D/v. Furthermore, by 
using d£/dt = mv dv / dt it follows from Eq. (IT^ that 
dvjdt ~ vva/{4:D). Inserting this (and ti, = 2£)/u) into 
Eq. (fT^ we obtain 

d f d 

-k(?^^/) =-4n^/diloss ■ (16) 

With the expression in Eq. (HI for i?ioss('c) a steady state 
solution {df /dt = 0) is then 


^ 1 f 2B{v-- 

/model OC o 6Xp 1 


Bn 


VVA 


(17) 


applicable for energies above e4)|| (or v > v^). 



FIG. 10: Regions in velocity space lost or gained by the 
velocity shell, wi < |v| < V2. 


E. Comparison of model for superthermal 
electrons to kinetic simulation data 

To compare Eq. (HZl to the kinetic simulation we use 
the energy distribution /e, here defined as a function of 
the relativistic kinetic energy £ = moc^i'y ~ !)• Using 
He = f fEd{-j — 1) = f fd?v, at non-relativistic energies 
we have oc (7 — 1)^/^/. Fig. ITT] shows Je for the 
points marked 1 - 3 in Fig. [51 with the most energetic 
distribution corresponding to point 3, the furthest point 
away from the x-line. For 77190^(7 — 1) < e$|| we approx¬ 
imately have /b oc (7 — 1 )^/^ corresponding to / ~ const 
for the near flat-top part of the distributions. For ener¬ 
gies above e4)|| the form in Eq. (113 represents a good 
approximation to the simulation data. 



FIG. 11: Electron energy distribution Je for the points 

marked PI - P3 in Fig. (6] For E < e4>|| the flattop part 
of the distribution is observed with /_b oc (7 — The 

dashed green/black line is /model in Eq. (I17II applicable for 
E > e4>||, evaluated with va/c = 0.08, e4>||/mc^ = 0.22 and 
B/Boo = 0.3. 

The distributions in Fig. [TT] demonstrate that signifi¬ 
cant electron energization is possible in an open exhaust 
of a single X-line reconnection configuration. We find 
that not only does $|| provide the electrons with an ini¬ 
tial energy boost as they enter the reconnection region. 
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the structure of also reduces the free-streaming along 
field lines in the reconnection exhaust. This permits the 
accumulation of energetic electrons, heated mainly by 
perpendicular electric fields during their repeated 

bounce-motion across the exhaust. Energization thereby 
becomes effective throughout the reconnection exhaust 
at much larger scales than the kinetic length scales of 
the electrons and ions. 


V. COMPARISONS OF NUMERICAL RESULTS 
WITH SPACECRAFT OBSERVATIONS 

A. Electron bulk energization during magnetotail 
reconnection 

Several decades of in situ spacecraft observations show 
that strong kinetic effects are present during magnetic 
reconnection in the Earth’s magnetosphere [13, Ha, 113 ■ 
Particularly relevant to the role of $|j are the electron dis¬ 
tributions documented by Nagai et al, which, similarly 
to Fig.|9lja,b), reveal the presence of cold beams directed 
towards the X-line while energized electrons move away 
from the reconnection region (see Fig. 4 in Ref. [13] )• The 
importance of these electron beams is also emphasized 
by measurements taken by the Geotail and the Cluster 
spacecraft [H, [13, Hll . Within the reconnection outflow 
the electrons often have a characteristic isotropic flat- 
top distribution, where the phase space density, /, of the 
electrons is nearly constant from thermal energies (tens 
to hundreds of eV) up to several keV. As detailed by 
Egedal et al. [l^, the distributions of Fig. IH] resemble 
closely distributions observed during reconnection events 
in the magnetotial. Furthermore, the extensive study of 
flat-top distributions by Asano et al. provides the spatial 
structure of the magnetotail exhausts where these flat- 
top distributions are observed (see Fig. 22 in Ref. [llj), 
which resemble the structure of <i>|| in Fig. [6]). 

The kinetic simulation results, the scaling law for $|| 
in Eq. m, and threshold in Eq. o for non-adiabatic 
electron behaviour all suggest that Peoo is an important 
parameter for the dynamics of the electrons during recon¬ 
nection. To explore the validity of these results to recon¬ 
nection in the magnetotail we study the inflow values of 
Pe as well as $|| inferred from a number of reconnection 
events encountered by the Cluster mission. 

The four Cluster satellites are moving together in coor¬ 
dinated polar orbits around the Earth, their internal sep¬ 
aration changing over the years of operation (year 2000 to 
present day). The data analyzed here is obtained through 
the Cluster Active Archive. The magnetic field data is 
provided by the Flux-Gate Magnetometer (FGM) exper¬ 
iment [4l| and the ion plasma data by the Cluster Ion 
Spectrometry (CIS) experiment [d^j. The electron data 
are from the Plasma Electron and Current Experiment 
(PEACE) [ 43 . In the analysis, only the PEACE data 
points with values above the background electron flux 
count that were also flagged with quality number 3 or 4 


(good for publication) were used. Both the HEEA (high 
energy sensor) and the LEEA (low energy sensor) data 
were included. To avoid photo electron measurements 
we consider only electrons at energy levels higher than 
70 eV. 

The data selected for our analysis were recorded dur¬ 
ing encounters with 21 separate magnetic reconnection 
events in the Earth’s magnetotail at about 18 Re. The re¬ 
connection events are listed in Ref. [13 along with the sig¬ 
natures of magnetic reconnection, as observed by space¬ 
craft flying through the reconnection region on a path 
parallel to the magnetotail neutral sheet. During all of 
the encounters, some or all of the spacecraft observed 
both reconnection outflow regions. The inflow regions 
were observed during some of the encounters. The signa¬ 
tures of both inflow and outflow regions are described in 
Ref. [13. It must be noted that while the outflow region 
is easily identified in the data, the inflow region is less 
distinct and this characteristic may introduce errors. 





FIG. 12: Characteristic electron distribution functions 
recorded by the Cluster Mission during a reconnection event 
on August 21, 2002. 

To determine the upstream electron characteristics for 
the observed magnetic reconnection encounters we con¬ 
sidered data recoded by spacecraft during excursions 
from the outflow regions into areas of lower plasma den¬ 
sity and stronger magnetic field. Still, these excur¬ 
sions did not bring the spacecraft into the magnetotail 
lobe. As an example of the data applied in this analysis, 
Fig. [m shows electron phase-space distributions recorded 
by Cluster 4 during a reconnection event encountered on 
August 21, 2002. The distribution in Fig. IT^ al is rep¬ 
resentative of the upstream reconnection region, and it’s 
thermal component is well described by a Maxwellian fit 
yielding inflow parameters of Tg ~ 80 eV and rie — 0.31 
cm“^. At the magnetic field of R ~ 30 nT, the nor¬ 
malized electron pressure for the inflow can then be esti- 
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mated as /3eoo — 0.011. 

The distribution in Fig. [T^ b'l is also for a location in 
the inflow region where the observed electron anisotropy 
develops in agreement with the kinetic model in Refs. |13l . 
fisj] . Meanwhile, the distributions in Figs. IT^ c.dl are ob¬ 
served in the exhaust and resemble the numerical distri¬ 
butions of Figs. in{a,b). Given this resemblance, a local 
value of e$|| can be identified as the energy of the incom¬ 
ing electron beams and/or the energy at which the flattop 
part of the distributions terminate. By considering all 
the measured distributions available, the maximal value 
of <i)|| for this event was estimated to be <i>|| ~ 9 keV, 
corresponding to a normalized potential e<i>||/Te ~ 110 
similar to the numerical values obtained above. 

To obtain the lobe plasma beta values we applied the 
electron and magnetic field data for the observations of 
the magnetotail lobe closest in time to the reconnection 
outflow observations. In table U we provide the values 
of $ 11 , j3\obe as well as the inflow values of rie, Te, B, 
and /3inflow(= Poo) estimated for the reconnection events. 
Out of the 21 events analyzed, we obtained inflow and 
lobe plasma beta values for 18 of them using data from 
the FGM and the PEACE instruments. Compared to the 
inflow values inflow characterizing the plasma feeding the 
reconnection region, the lobe values /3iobe were typically 
lower by an order of magnitude. 

The results from the described analysis of the Cluster 
data are summarized in Fig. 1131 which provides the values 
of e^ii/Te as a function of Peoo- For comparison, the 
red symbols represent the range of e$||/Te recorded in 
the numerical simulations (with values from particular 
figures above indicated). In addition, the blue line is 
e$||/Te of Eq. ([^I, which yields a lower bound for the 
observations and simulations. 

The range of e$||/Te from the simulation presented 
here (with Peoo = 0.003) is in good agreement with the 
spacecraft observations. However, the spacecraft data 
suggest that the transition to the regime influenced by 
electron holes and double layers (large values of e$||/Te) 
occurs at Peoo — 0.1. In contrast, all simulations we have 
studied with — 0.1 are well described by Eq. © 
(the blue line) [^. Nevertheless, the threshold for the 
transition to non-adiabatic inflow electron, P^oo < 0.02 
in Eq. 0, is in reasonable agreement with the spacecraft 
observations. 


B. Superthermal electrons in the magnetotail 

In a previous study of magnetotail observations we 
found that superthermal electrons often acquire a con¬ 
stant energy-gain, A£, independent of their initial energy 
(45l| . As an example, in Fig. 1141 we provide measurements 
of superthermal electrons recorded by the Cluster Mis¬ 
sion durin g th e much studied October 1, 2001 reconnec¬ 
tion event [^. Time series measurements by the RAPID 
instrument are given in a). For the selected time points 
for which the distributions are shown in b), we apply a 



FIG. 13: Values of e’Fn/Te as afunction of deoo- The black cir¬ 
cles are calculated based on the values in Table[T]inferred from 
Cluster spacecraft data. The red symbols represent values of 
kinetic simulations. The values of d’y/Te for deoo = 0.05 and 
0.15 are from simulations presented in Ref. |^. The blue line 
is the adiabatic model in Eq. ©• This model does not include 
the effects on 4>|| caused by electron holes and double layers, 
and thus represents a lower theoretical bound for e4>||/re. 


TABLE I: Plasma beta and guide field values for magnetic 
reconnection inflow regions and the lobe. 


Date 


$11 /kV Aobe 

Te/eV 

Ue/cm ® 

B/nT Anflow 

2001 08 

22 

8 

0.003 

210 

0.07 

26 

0.008 

2001 09 

10 

0.8 

0.003 

340 

0.36 

15 

0.21 

2001 09 

12 

5 

0.001 

90 

0.25 

20 

0.026 

2001 10 

01 

10 

0.008 

150 

0.28 

30 

0.018 

2001 10 

08 

2 

0.004 

500 

0.45 

20 

0.22 

2001 10 

11 

2 

0.03 

650 

0.14 

15 

0.15 

2002 08 

21 

9 

0.003 

80 

0.32 

30 

0.011 

2002 08 

28 

2.1 

0.0003 

200 

0.2 

20 

0.038 

2002 09 

13 

3 

0.001 

220 

0.23 

24 

0.034 

2002 09 

18 

11 

0.0009 

230 

0.24 

20 

0.054 

2002 10 

02 

11 

0.002 

170 

0.064 

12 

0.030 

2003 08 

17 

5 
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Liouville mapping technique to obtain the energy gains 
shown in c). Again, more details of this analysis are giv¬ 
ing in Ref. m. 

For all the events analyzed it is found that, for the 
most energitic electrons, AS is constant, independent of 
the initial energy. In Ref. [1^ we then concluded that the 
energization of the superthermal electrons is set by e$||, 
but our understanding of this has changed and is now 
very different. From the above analysis it is clear that it 
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is the bulk-energization (including flat-top distributions) 
that is controlled by e<i>||, and the observed superthermal 
electrons energization is better described by Eq. (|T^ . al¬ 
lowing for energization to energies much larger than e$||. 
However, because of the initial hard power-law spectrum 
f{v) oc v~^ of incoming lobe electrons and because A£ in 
Eq. (|12l) is proportional to the initial energy, the present 
2D model predicts superthermal electrons in excess of 
the observations. This discrepancy can be explained by 
the finite y-extent of the magnetotail. In a realistic ap¬ 
plication, the effectiveness of the present model needs to 
be limited, as the maximum energization level is set by 
the potential drop along the extent of the reconnection 
X-line. The constant levels of A£ deduced in the analy¬ 
sis of Ref. (dHj i are then consistent with our superthermal 
heating mechanism in combination with a finite dusk- 
dawn (or for observations: dusk - spacecraft) potential 
drop during magnetotial reconnection. 


Cluster 3, Time in Sec since 2001-10-01, 09:48:00 UT 




logl0(E/eV) 



FIG. 14: Electron phase-space densities from the RAPID 
measurement on Cluster 3. (b) Electron distribution observed 
at separate time points, (c) Spectra of A£ for three selected 
times. The observed spectra of A£, independent of £ for large 
£, are consistent with Eq. m when taking into account the 
finite extent of the systems in the magnetotail y-direction. 


C. Solar flares 

Compelling observational evidence exists for confine¬ 
ment of energetic electrons by parallel electric fields dur¬ 
ing solar flare events. This evidence is summarized in pa¬ 
pers by Li. et al, [H, who (as mentioned above) 

explored the role of double layers for reducing the free 
streaming losses of reconnection energized electrons along 
the field lines. The evidence includes the formation of 
energetic electron populations detected in the vicinity of 
the expected loop-top reconnection sites. The lifetime 


of these energetic populations are inferred to be much 
longer than the thermal electron transit time a , sug¬ 
gesting that parallel electric fields must be important to 
reduce free streaming losses along magnetic field lines. 



FIG. 15: Schematic model of a solar flare 

Based on the Masuda flare model [i^. Fig. [15] provides 
a schematic illustration of the electric and magnetic ge¬ 
ometry suggested by the observations and required for 
confinement of electrons by parallel electric fields. Again, 
the analysis by Li. et al. and also earlier authors (see ref¬ 
erences in Ref. 0 ) have already suggested that electric 
fields in double layers may help confine the energized elec¬ 
trons. In addition, the present analysis suggests that the 
parallel electric fields are not only important for confining 
the electrons but are an integral component to the overall 
energization processes. For example, the observed and re¬ 
markably efficient energization of the bulk electrons from 
about IkeV up to lO-lOOkeV [I| is consistent with the 
development of an acceleration potential with a magni¬ 
tude of e$||/Te = 100 seen in the magnetotail and in 
our simulation. In addition, the most energetic electrons 
reaching the MeV energy range can also be accounted for 
by Eq. (HZl). In fact, for energies just above e<i>|| , Eq. (El) 
yields a hard powerlow spectrum f(v) oc u ^ for about 
an order of magnitude in energy, sufficient to explain the 
most energetic electrons in the solar flare observations. 
We also note that the /3eoo-threshold {f3eoo < 0.02) for 
the non-adiabatic electron dynamics to occur is easily 
satisfied in solar flare events 113. 
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In our model, the electrons are only energized by $|| 
as they enter the reconnection exhaust. Meanwhile, ions 
will be accelerated away from the reconnection region by 
these parallel electric fields. The path of the ions will 
not coincide with instantaneous magnetic field lines, and 
a fraction of the ions are likely to endure acceleration by 
to much larger values than e<I>|| during their transit 
through the reconnection exhaust. This may account for 
the energization of ions in the range of 10 MeV recorded 
by the RHESSI spacecraft [i^. Of course, it must be 
noted that our numerical simulation covers a domain size 
of about 300di, whereas the size of a solar flare is about 
lO^di. Furthermore, solar flares are believed to include a 
guide magnetic field of order unity, whereas the present 
simulation is for anti-parallel reconnection. Nevertheless, 
we expect that the described general processes leading to 
non-adiabatic parallel electron dynamics and strong elec¬ 
tron energization will also be applicable to reconnection 
in very large systems. 


VI. CONCLUSIONS 

The heating by E-v^ has been found generic to plasma 
flows driven by magnetic tension and is common to 
a range of models for electron energization. An example 
is Drake’s Fermi acceleration model considering a bath 
of reconnecting magnetic islands . Likewise the term 
E • V£) is also responsible for driving powerlaw distri¬ 
butions in highly relativistic pair plasmas (50j |. A third 
example is the energization documented while propagat¬ 
ing test-particles through the fields of Ideal-MHD simula¬ 
tions [5ll |. where confinement may occur when electrons 
bounce between dipolarization fronts and the stronger 
magnetic fields close to Earth [Hill. However, Ideal- 
MHD assumes Ey = 0 and the described effects of local¬ 
ized electron trapping and the initial energization by Ey 
are therefore omitted. 

The key new feature presented here is the importance 
of the magnetic field aligned electric fields Ey for initial 
energization of the electrons. Because the subsequent 
energization by E_l is propotional to £, the initial en¬ 
ergy boost that electrons acquire from e$|| largely deter¬ 
mines the overall efficiency of the energization process. 
The large amplitude structure of <I>y also limits the free 
streaming particle losses along magnetic field lines, which 
is essential for confining the electrons and shaping the 
spectra of their energy distribution. The heating is effec¬ 
tive in the full exhaust of a single X-line configuration, 
and does not require pressure anisotropy to develop. 

The large values of Ey are observed within density cav¬ 
ities and develop in association with the strong double 
layers. A threshold, /3eoo < 0.02 is derived for the strong 
double layer formation, suggesting that these are likely to 
be present during solar flare events. Indeed, the level of 
electron energization observed in the simulation is consis¬ 
tent with that observed during solar flares. Furthermore, 
we show the details of the electron distribution functions 


described in our analysis are in agreement with in situ 
spacecraft measurements obtained over the last decade 
during reconnection events in the Earth’s magnetotail. 
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Appendix A: Use of the guiding center model for 
estimating the electron energization rate 


In this appendix we first discuss why the guiding center 
model can provide an accurate prediction of the electron 
heating rate for the pitch angle mixed exhaust. We then 
estimate the heating rate based on this model and recover 
result given above in Eq. (HID- We thereby validate the 
above estimate and highlight the similarity (and differ¬ 
ences) between the results in the present paper compared 
to those of earlier works using the guiding center model 


From Fig. [T^b) above it is clear that the main heating 
source of the electrons is the perpendicular electric fields, 
and thus, the local heating results from the term eneVe_L- 
E, where Vej_ is the electron fluid velocity also appearing 
in the electron momentum equation: 


E-fVe xB = -V-Pe/(ene) . (18) 


Due to the low values of B in the exhaust, the mag¬ 
netic moments of the electrons are not conserved. The 
particle motion therefore becomes chaotic and all pitch 
angle information is lost as electrons pass through the 
midplane. We are therefore only interested in the av¬ 
erage electron behavior as a function of energy. It is 
clear that the guiding center approximation with the 
drifts vd = Vfe -f vvs does not formally apply (where 
Vfe = me'CyB X K/{eB^) with k = b • Vb is the cur¬ 
vature drift and vys = me^^B x VE/(2eE^) is the 
gradient-B drift). Nevertheless, by direct evaluation of 
the expression eneVej_ = e f v/jfd^v + VxM, with the 
magnetization eM = b J /if d^v, it is well known that 
V e_L in Eq. (flSl) is recovered for the relevant limit where 
Pe has no off diagonal stress [U- Therefore, although 
the guiding center approximation does not account for 
the chaotic motion of the individual electrons, it does ac¬ 
curately predict the total fluid drift. Furthermore, the 
predicted contribution to Ve_L from particles within a 
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given energy interval is also accurate. The guiding cen¬ 
ter approximation thus allows us to estimate the average 
electron drift (and associated energization) as a function 
of energy. 

To derive an expression for d£/dt using the guiding 
center model, we note that previous studies in similar ge¬ 
ometries have found that curvature drifts dominate the 
particle motion in the direction of the reconnection elec¬ 
tric field [S^. Therefore, we here estimate the energy 
changes over a single bounce orbit caused by the cur¬ 
vature drift, AS = e /J"’’• E^ecdt. Using |«;| = 1/Rc 
where Rc is the radius of curvature for the magnetic field 
line and T}, ~ ttRc/v^ we obtain AS ~ ifrecTTWelun |/i? ~ 
Ej-ecinv/B, as the average energization per bounce. Us¬ 


ing i?rec — waT^o/IO we obtain 

dS AS vaBq mv va 9 

— = -~-~ - mv , 

dt n 10 Bn lorb 

where we have used B ~ Bq/IQ. Thus, we find that 


dS 

dt 



(19) 


This order of magnitude estimate of the heating rate is 
consistent with the result in Sec. IIII PI In fact, taking 
the orbit length as four times the half exhaust width, 
^orb = 4Z1, we recover the expression in Eq. 
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